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Electron  Paramagnetic  Resonance  of 
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L.  Rimai  and  G.  A.  deMars 

Raytheon  Company 
Research  Division 
Waltham,  Massachusetts 

Abstract 

This  paper  reports  electron  paramagnetic  resonance  data  of  substitutional 

Gd  impurities  in  the  cubic  and  tetragonal  phases  of  barium  and  strontium 

titanates.  The  pertinent  spin  Hamiltonian  parameters  are  determined  as  a 

function  of  temperature  in  these  two  mentioned  phases  of  both  crystals.  The 

temperature  for  the  cubic  to  tetragonal  phase  transition  in  SrTi03  is  determined 

from  these  data  to  be  (110±2.  5°K);  a  result  that  has  been  checked  by  the 

observation  of  the  temperature  dependence  of  the  birefringence  pattern  of 

pure  (undoped)  crystals.  This  experiment  is  also  briefly  described.  In  the 

latter  crystal  we  find  a  smooth  transition,  contrary^,  to  the  corresponding 

+  +  + 

case  in  BaTi03.  The  Gd  resonance  results  are  compared  with  previous 
data  on  EPR  of  iron  group  impurites,  and  some  comments  and  conjectures 
are  offered  as  to  possible  correlation'  of.these  results  with  the  temperature 
dependence  of  certain  properties  of  the  host  crystals  in  the  neighborhood  of  the 
phase  transition. 


I.  INTRODUCTION 

The  purpose  of  this  paper  is  to  present  and  discuss  results  obtained 
by  examining  the  electron  paramagnetic  resonance  (EPR)  of  Gd  occurring 
as  a  substitutional  impurity  in  the  barium  and  strontium  titanates.  A  number 
of  publications 1  ~ ^  have  dealt  with  the  EPR  of  iron  group  impurities  in  these 
crystals.  Size  considerations  indicate  that  such  ions  would  preferentially 

XJ.11 

substitute  at  the  Ti  site  whereas  larger  ions,  e.  g.  ,  the  rare  earths, 

would  probably  substitute  for  the  alkaline  earth.  We  decided,  therefore, 

7 

that  an  investigation  of  a  rare  earth  impurity  would  be  useful.  Among 

+++  ++ 

the  rare  earths,  only  the  S  state  ions  Gd  and  Eu  are  susceptible  to 
investigation  over  a  wide  range  of  temperatures. 

g 

It  is  well  established  that  BaTiO^  exhibits  an  abrupt  phase  transition 

in  the  neighborhood  of  120°C  from  the  high  temperature  cubic,  unpolarized 

phase  to  a  tetragonal,  ferroelectric  phase  which  is  the  stable  phase  at  room 

temperature.  This  phase  transition  shows  its  effect  quite  obviously  on  the 

+ ++  2 

temperature  dependence  of  the  EPR  of  the  Fe  impurities.  Our  results 

111 

show  a  similar  effect  in  the  EPR  of  Gd  .  A  sharp  change  is  observed  at 
the  transition  temperature. 

For  a  long  time,  SrTiO^  was  considered  to  be  cubic  at  all  temperatures. 
EPR  studies  of  Mn  impurities^  indicated,  however,  the  presence  of  a  phase 
transition  in  the  neighborhood  of  the  temperature  of  liquid  air.  One  cannot 
be  sure  in  this  case  that  the  Mn  v  ion  is  perfectly  substitutional  in  either 
phase.  This  makes  the  analysis  of  these  results  insofar  as  the  characteriza¬ 
tion  of  a  phase  transition  in  the  host  somewhat  difficult.  The  most  reasonable 
interpretation  indicates,  however,  that  there  is  a  structural  change  from  a 
cubic  phase  stable  at  the  higher  temperatures  to  a  tetragonal  phase  stable 
below  the  boiling  point  of  N2- 

++  + 

EPR  studies  of  Fe  in  SrTiO^  show  also  that  it  is  a  substitutional 

impurity  with  cubic  site  symmetry  at  room  temperature,  but  with  a  small 

4  5 

tetragonal  distortion  present  at  77°K  and  below.  ’  A  more  detailed  tempera¬ 
ture  dependence  of  the  spectrum  in  this  case  could  presumably  have  character- 
i£®. d  the  transition  in  more  detail. 

Evidence  for  a  low  temperature  phase  transition  in  SrTiC^  from  dielectric 
studies  'has  been  presented  earlier.  This,  however,  gave  rise  to  some 
controversy.  See  references  to  work  by  Graincher  as  mentioned  in  refs,  4  and  8 
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The  present  results,  in  particular  the  detailed  temperature  dependence 

of  the  Gd  spectrum  in  SrTiO^,  indicate  the  presence  of  a  cubic  to  tetragonal 

phase  transition  close  to  110°K.  Within  the  limits  of  experimental  accuracy, 

this  is  a  smooth  transition,  in  contrast  to  the  sharp  transition  found  in  BaTiO^. 

In  order  to  obtain  definite  evidence  that  this  transition  is  not  a  property  of 

the  neighborhood  of  the  impurity  alone,  we  have  examined  the  temperature 

dependence  of  the  birefringence  pattern  of  single  crystals  of  pure  SrTiO^  and 

determined  that  a  phase  transition  in  these  crystals  does  indeed  occur  at  this 

temperature.  More  recently  additional  evidence  for  this  transition  has  been 
9  10 

obtained.  "  Actually,  from  the  fact  that  throughout  this  transition  the  EPR 
lines  remain  narrow  is  in  itself  evidence  that  this  transition  is  brought  about 
by  a  cooperative  phenomenon  in  the  host  crystal.- 

Some  comment  on  the  birefringence  experiments  will  be  given  in  Section  II, 
together  with  a  discussion  of  the  experimental  techniques  involved  in  the  EPR 
spectroscopy.  The  results  of  the  EPR  measurements  will  be  presented  in 
detail  in  Section  in.  Section  IV  is  devoted  to  a  discussion  of  our  results  in 
light  of  present  theories  of  ferroelectrics  of  the  BaTi03  type. 

We  feel  that  these  comparative  studies  of  the  tetragonal -cubic  transition 
in  the  two  crystals  considered  should  be  of  help  in  bringing  about  a  more  com¬ 
plete  understanding  of  the  properties  of  either. 


H.-  EXPERIMENTAL  TECHNIQUES 

A.  Samples  and  Crystallographic  Considerations 

The  SrTi03  samples  were  purchased  from  the  National  Lead  Company 

and  contained  about  0.  l°/o  atomic  Gd.  The  BaTiO^  crystals  were  grown  in 

this  laboratory  by  Mr.  B.  diBenedetto  using  flux  techniques.  Crystals  contain- 

ing  atomic  concentrations cf  Gd  from  0.2  to  0.01  /o  were  examined.  The 

less  concentrated  crystals  were  slightly  colored  plates  with  their  larger  faces 

along  (100)  planes,  although  (110)  faces  were  also  present  on  some.  The  crystals 

+++ 

with  higher  Gd  concentrations  exhibited  generally  the  same  geometry  as  the 
others;  however,  they  contained  in  most  cases  dark  occlusions  of  unreacted  flux 
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material  and  were  appreciably  lossier.  Nevertheless,  the  same  spectrum  was 
encountered  in  all  the  samples  examined  aside,  of  course,  from  intensity  varia¬ 
tions  due  to  concentration  differences.  Our  temperature  resolution  (±1-  5°C)  in 
the  region  above  room  temperature  whs  not  sufficient  to  measure  the  slight 
impurity  concentration  dependence  of’the  transition  temperature  which  may  be 
observed  by  other  means. *  * 

In  the  high  temperature  phase,  both  crystals  have  the  cubic  perovskite 
8 

structure,  where  there  is  only  one  metal  ion  of  each  kind  per  unit  cell.  Thus 

+++ 

the  metallic  sites  have  cubic  point  symmetry,  and  for  substitutional  Gd  the 
spectrum  exhibits  this  symmetry  with  axes  parallel  to  edges  of  the  unit  cell. 

A  number  of  subsidiary  weak  lines  are  observed  in  the  spectrum  of  Gd  in 
SrTiO^.  These  lines  may  be  attributed  to  a  small  fraction  of  Gd  ions  which 
do  not  enter  as  perfect  substitutional  ions,  conceivably  due  to  some  axial  symmetry 

7 

in  the  mechanism  of  charge  compensation  for  the  valance  unbalance.  We  shall 
not  be  concerned  with  such  spectra  in  the  present  paper. 

In  the  tetragonal  phase,  twinning  of  the  crystal  takes  place  resulting  in 
the  formation  of  domains.  There  are  three  equivalent  types  of  domains 
according  to  which  of  the  three  [  100]  directions  is  parellel  to  the  tetragonal 
distortion.  Thus,  with  an  arbitrary  orientation  of  the  magnetic  field  H  with 
respect  to  the  crystal  axes,  three  different  spectra  should  be  observed. 

However  if  H  is  confined  to  a  [  110]  plane  this  number  is  reduced  to  two. 

Due  to  the  simplification  thereby  achieved,  that  is,  the  manner  in  which  the 

/ 

angular  dependence  data  on  SrTiO^  were  taken. 

Since  the  BaTiO^  exhibited  natural  [  100]  faces  these  planes  were  used 
to  obtain  rotation  patterns.  In  this  case  there  was  a  strong  preferential 
orientation  of  the  domains  (their  tetragonal  distortion)  normal  to  a  particular 
face.  Thus  rotation  in  a  natural  [  100]  face  parallel  to  the  direction  of  preferential 
orientation  yielded  the  best  signal -to -noise  ratio.  Lines  due  to  the  other 
domains  were  observed  also. 

By  the  term  ’’domains"  we  simply  mean,;,  in  this  case,  regions  of  uniform 
tetragonal  distortion.  Their  boundaries  are  not  necessarily  determined  by  the 
long  range  dipolar  interaction.  They  are  more  probably  related  to  the  distribution 
of  strains  in  the  crystal. 
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B.  Remarks  on  EPR  Experiments 

The  EPR  data  were  taken  using  a  spectrometer  built  to  cover  the  Ku 

microwave  band  (12-18  kMc/sec).  The  system  was  of  fairly  standard  design 

using  a  balanced  bolometer  detection  scheme  with  25  c/ sec  field  modulation 

12 

phase  sensitive  narrow  banding.  The  klystron  was  frequency  stabilized  to 
the  sample  cavity  resonance.  This  was  quite  necessary  in  view  of  the  high 
and  appreciably  temperature  dependent  dielectric  constant  of  the  samples. 

In  the  temperature  range  from  300°  to  77°K  a  platinum  thermometer  imbedded 
in  the  cavity  wall  was  used.  The  accuracy  in  sample  temperature  thereby 
achieved  was  ±  0.  5°K.  '  In  the  range  above  300°C  the  temperature  was  measured 
by  means  of  a  thermocouple. 

An  interesting  effect  appears  in  the  neighborhood  of  the  transition  for 
BaTiOy  When  the  transition  temperature  is  approached  from  below  the  tetragonal 
spectrum  disappears  and  the  cubic  spectrum  appears.  There  is  a  small 
hysteresis  region  in  which  the  two  spectra  should  exist  simultaneously.  This 
temperature  interval  is  very  small-(about  2°C.  )and  we  have  not  observed  it 
due  to  lack  of  resolution.  We  observed:,  however,  that  as  the  transition  was 
approached  there  was  a  marked  loss  in  the  intensity  of  the  tetragonal  spectrum 
and  also  the  cubic  spectrum  when  it  first  appeared  was  quite  weak  increasing 
in  intensity  with  T  right  above  the  transition.  There  was  definitely  a  large 
loss  of  total  intensity  near  the  transition  temperature.  There  was^  howeverr-no 
observable  increase  of  the  losses  in  the  sample  as  shown  by  the  fact  that  the 
coupling  and  Q  of  the  cavity  did  not  change.  This  effect  seems  to  indicate  a 
strong  increase  in  the  line  widths  of  all  but  the  central  (l/ 2,  -  l/ 2)  EPR  lines 
when  the  transition  temperature  is  approached.  This  in  turn  might  be  a  con¬ 
sequence  of  ionic  motion  in  the  sample  of  frequency  not  too  far  removed  from 
that  corresponding  to  the  tetragonal  microwave  splitting. 

C.  The  Birefringence  Experiment 

When  viewed  between  crossed  polaroids  pure  SrTiO^  crystals  exhibit  a 

complex  birefringenece  pattern  even  at  room  temperature.  When  a  plate  (about 

1  mm  thickness)  is  examined  with  monochromatic  light  it  exhibits  distinct 

regions,  each  with  the  characteristic  45°  angular  separation  between  directions 

*  The  accuracy  with  which  we  are  able  to  state  the  temperature  for  the  phase 
transition  in  SrTi03  is  less,  however.  This  is  due  to  the  limitation  imposed  on 
spectral  resolution  by  the  width  of  the  absorption  lines. 
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of  light  polarization  for  maximum  and  minimum  transmission.  These  regions 
are  separated  by  lines  that  are  always  opaque  (non -birefringent) .  Such  a  pattern 
may  be  explained  by  the  existance  of  uniformly  strained  regions  in  the  initially 
cubic  crystal.  All  crystals  examined  exhibited  such  patterns  and  annealing 
for  a  few  hours  at  1600°C  was  not  effective  in  altering  them. 

These  strain  patterns  were  examined  as  a  function  of  the  temperature  between 
300°K  and  77°K  with  the  crystal  mounted  in  an  optical  dewar.  No  noticeable 
changes  occurred  in  the  pattern  between  300°  and  130°K.  Below  100°K  a  tempera¬ 
ture  independent  pattern  was  observed  again,  markedly  different,  however,  from 
the  room  temperature  pattern.  As  the  sample  was  slowly  warmed  up  from  100°K 
the  dark  boundaries  between  the  birefringent  regions  started  to  move  and  broaden 
until,  in  the  neighborhood  of  110°K,  the  crystal  was  opaque.  At  this  temperature 
rotation  of  the  polaroids  showed  the  crystal  to  be  essentially  optically  isotropic. 

As  the  temperature  continued  to  rise  new  transparent  regions  started  to  appear 
and  grow  until  around  130°K  a  pattern  appeared  which  was  identical  to  what  was 
observed  at  room  temperature. 

Even  if  the  birefringence  pattern  in  the  crystal  above  and  below  110°  is 
ascribed  to  strains,  the  rearrangement  in  this  pattern  that  occurs  at  this  tempera¬ 
ture  definitely  indicates  a  structural  modification,  in  the  crystal.  This  result 
gives  us  confidence  in  interpreting  the  transition  of  the  EPR  spectrum  as  a 
result  of  a  structural  change  in  the  host  crystal.  Thus  the  information  obtained 
from  the  EPR  may  be  now  used  as  a  description  in  some  detail  of  the  nature  of 
this  transition. 


III.  EPR  RESULTS 

7  3  |  A 

The  ground  state  of  the  free  Gd  ion  is  4f  S,  To  describe  the  energy 
levels  resulting  from  the  splitting  of  this  degenerate  state  by  the  combined  effect 
of  the  crystalline  electrostatic  field  and  the  external  magnetic  field  H,  we  write 
a  spin  Hamiltonian^ 

3C(S)  =  gpH-S+X  (-Dmb  Y  m(S) 

!,m  £’m 


(1) 
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S  is  the  spin  vector;  S  =  7/2  of  the  free  ion  ground  state.  The  first  term  is 

the  Zeeman  energy  with  an  isotropic  splitting  factor  g.  Crystalline  anisotropy 

14 

of  g  is  usually  negligible  for  S  ground  states.  The  second  term  represents 
the  interaction  of  the  spin  S  with  the  electrostatic  field  at  the  site  that  the  ion. 
occupies  in  the  crystal  unit  cell.  We  expand  this  interaction  in  the  spin  dependent 
operators  Y.  (S)  which  have  the  property  that  under  rotation  they  transform 

a  f  in  y 

as  the  X,mth  spherical  harmonic.  Using  such  an  expansion, the  matrix 

elements  of  the  crystalline  field  Hamiltonian  between  eigenstates  of  S  =  Mq 

|  'y  z  o 

may  be  readily  obtained  from  tables  of  3j  symbols,  since  the  reduced  matrix 

elements  which  are  in  this  case  diagonal  with  S  =  7/2  may  be  incorporated  in 

the  constants  b  .  Furthermore,  we  normalize  the  Y,  (S)  in  such  a  way 
jl  f  m  jl  t  m 

that  the  diagonal  matrix  elements  become  the  smallest  possible  integers.  This 

is  achieved  by  incorporating  into  the  b  factors  depending  only  on  i.  Using 

this  recipe  to  write  the  matrix  elements,  the  coefficients  b  that  we  deter- 

mine  from  our  data  are  in  the  same  convention  as  those  for  which  the  same 

14  18  19 

notation  has  been  generally  used.  ’  ’ 


The  expression  (1)  applies  in  a  general  coordinate  system.  For  a  particular 
point  symmetry  we  may  choose  axes  for  which  the  largest  number  of  coefficients 
b.  vanish.  If  we  have  tetragonal  symmetry,  and  if  we  choose  axes  along  the 

X  f  IH 

three  [  100]  directions  with  z  the  fourfold  axis, 


X  -  gpH*  S  +  b2joY2,0  +  b4,0Y4,0  +  b4,4(Y4,4  +  Y4, -4)  +  b6,0Y6,0 


+  b6,4(Y6,4  +  Y6,-4} 


(2) 


We  have  restricted  (2)  to  terms  with  &  £  6,  since  higher  i  terms  reduce 
to  linear  combinations  of  the  above  for  S  =  7/2.  In  the  particular  case  of  cubic 
symmetry  the  number  of  independent  coefficients  is  decreased  since  the  follow¬ 
ing  relations  hold:* 


0 


b 

14  °4, 0 


\/l4  , 

6,4  ~  "  2  °6, 0 


(3) 
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If  the  magnetic  field  is  parallel  to  z,  and  if  we  choose  a  representation 

for  which  S  is  diagonal,  (2)  will  be  diagonal  except  for  the  M  =  4  terms.  It 
z 

turns  out  fortunately,  however,  that  the  frequency  at  which  these  experiments 
were  performed  (w/2 ir  =  16  kMc/sec)  is  sufficiently  high  so  as  to  permit  the 
treatment  of  the  crystal  field  interaction  by  perturbation  theory;  first  order  for 
terms  in  i  =  4.  and  6,  second  order  at  the  lower  temperatures  for  l  =  2. 

Thus  from  data  obtained  with  H  parallel  to  the  tetragonal  axis  the  three 
constants  b^  q,  b^  q,  and  b^  ^  may  be  determined  from  first  order  formulae 
alone.  The  g  value  is  determined  from  the  position  of  the(-l/2  ,  l/2) transition 
also  from  these  data. 


For  a  tetragonal  spectrum,  b^  ^  and  b^  ^  are  obtained  from  data  taken 
when  H  is  in  a  plane  normal  to  the  tetragonal  axis,  either  parallel  to  a  (100) 
or  to  a  (1 10). 


Table  I  contains  a  collection  of  values  of  these  parameters  obtained  at  a 
few  selected  temperatures.  Using  these  parameters, we  have  calculated  by 
perturbation  theory  as  outlined  above  the  angular  dependence  of  some  of  the 
transitions  and  obtained  satisfactory  agreement  with  experimental  results  as 
shown  in  Fig.  1,  2,  and  3.  The  verification  of  the  angular  dependence  serves 
the  purpose  of  excluding  the  presence  of  extra  terms  in  the  spin  Hamiltonian, 
e.g.  Y2  2  +  due  to  a  rhombic  distortion  or  of  a  tetragonal  term,  Y.,q,_ 

in  the  cubic  phase  which  would  mean  that  the  impurity  is  not  substitutional. 

The  comparison  of  the  values  in.  Table  I  among  themselves  and  to  the 
Zeeman  energy  for  16  kMc/sec  justifies  the  use  of  perturbation  theory  as 
outlined.  As  a  consequence,  we  do  not  have  appreciable  mixing  of  different 
|Mg>  states  and  thus  relative  line  intensities  serve  to  identify  the  various 
transitions.  To  obtain  good  agreement  with  ratios  of  rf  transition  probabilities 
for  the  BaTiO^  data,  integrated  intensities  had  to  be  used  since  there  is  notice¬ 
able  broadening  of  these  lines.  This  broadening  increases  as  we  move  away 
from  the  central (l/2,  -l/2)line,  presumably  due  to  fluctuations  in  the  crystal 
field  parameters  throughout  the  imperfect  crystal. 

The  sign  of  the  parameters  for  BaTi03  has  only  relative  meaning.  However, 
for  SrTiO^  the  signs  are  absolute  since  the  fact  that  the  tetragonal  spectrum  was 
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observable  down  to  1. 5°K  permitted  the  direct  determination  of  the  sign  of  b^  q 
relative  to  which  we  have  taken  the  signs  of  the  other  parameters.  We  observed 
that  in  cooling  to  1.  5°K  the  intensity  of  the  outermost  high  field  (7/2  ~  5/2) 
transition  decreased  relative  to  the  corresponding  low  field  transition.  Thus 
b^  q  <  0.  The  assignment  of  absolute  signs  to  q,  b^  q  in  the  cubic  phase 
rests  on  the  assumption  that  b4  q  is  a  smooth  function  of  the  temperature  in 
the  neighborhood  of  the  transition.  This  is  reasonable  since  its  absolute  value 
is  certainly  so  and  a  reversal  in  sign  at  the  transition  without  change  in  absolute 
value  seems  unlikely. 

The  above  discussion  has  dealt  with  the  determination  of  the  values  of 
the  parameters  in  the  spin  Hamiltonian.  The  most  interesting  results  of  the 
present  work  are  the  temperature  dependence  of  the  EPR  spectrum  and  the 
temperature  dependence  of  the  parameters,  especially  that  of  b^  q. 

Figure  4  shows  the  temperature  dependence  of  the  observed  spectrum 
in  SrTiO^  when  H  is  parallel  to  the  tetragonal  axis  of  certain  domains. 

We  see  below  110'K  two  spectra  due  to  the  presence  of  ions  with  tetragonal 
axes  along  the  other  two  [  100]  directions.  We  also  observe  the  seven  transitions 
of  a  practically  temperature  independent  spectrum  above  110°K,  corresponding 
to  one  type  of  ion  only. 

Figure  5  shows  the  temperature  dependence  of  the  spectrum  in  BaTiO^ 
with  H  parallel  to  the  preferred  tetragonal  axis.  In  the  cubic  phase, due  to 
the  large  line  widths»the  various  transitions  could  not  be  resolved.  From  latent 
structure  in  the  lines, the  value  for  b4  q  shown  on  Table  I  is  deduced.  The  abrupt 
nature  of  the  transition  is  clearly  observable. 

The  parameter  with  the  strongest  temperature  dependence  is  b2  q.  This 
is  shown  in  Fig.  6  for  SrTiO^.  Possible  implications  of  this  result  will  be 
discussed  in  the  next  section.  Included  in  Fig.  6  is  the  temperature  dependence  of 

b  which  is  much  slower  than  that  of  b_  n. 

4,0  2,0 

The  values  of  b^  q  in  most  cases  are  very  small,  bordering  experimental 
error  limits,  and  will  not  be  further  discussed. 
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IV.  DISCUSSION 

Certain  aspects  of  the  above  results  are  of  interest  in  their  relation  to  the 
properties  of  the  host  in  the  vicinity  of  the  phase  transition.  In  what  follows  we 
discuss  some  of  these  points. 

A.  The  existence  of  the  cubic-tetragonal  phase  transition  in  BaTiO^  has 

been  established  and  (studied  by  many  techniques  previous  to  the  first  EPR 

investigations.  In  the  case  of  SrTiO^,  however,  the  initial  evidence  for  the 

existence  of  such  a  transition  in  the  neighborhood  of  1 1 0°K  has  been  obtained  from 
3  4 

EPR  results.  ’  In  spite  of  the  fact  that  in  such  experiments  one  is  dealing  with 
properties  of  impurity  ions  included  dilutely  in  the  host,  these  results  alone  permit 
us  to  conclude  the  existence  of  such  a  phase  change. 

If  one  had  assumed  the  crystal  to  be  cubic  at  all  temperatures,  the  tetragonal 
symmetry  of  the  spectra  at  low  temperatues  could  be  explained  by  saying  that  the 
impurity  is  not  perfectly  substitutional.  For  both  Gd+++  and  Fe+++,  a  reason  for 
this  c  ould  be  found  in  the  necessity  of  compensating  the  deviation  of  local  charge 
neutrality  due  to  the  difference  in  valence  between  these  ions  and  either  metallic 
component  of  the  host.  To  maintain  electrostatic  neutrality  the  impurity  ion 
would  move  away  from  the  substitutional  position  of  cubic  symmetry,  say  by  a 
small  displacement  parallel  to  the  [100]  direction.  Since  there  are  three  equi¬ 
valent  such  directions,  one  would  obtain  in  general  three  inequivalent  spectra. 

With  increasing  temperature  a  thermally  activated  motion  characterized  by  a  corre¬ 
lation  time  t  =  r  :exp(**A/kT)  would  enable  the  ion  to  jump  between  the  equivalent 
°  20 
positions  with  distortions  along  the  three  [  100]  directions.  Such  a  motion  would, 

upon  heating,  progressively  average  out  the  non- cubic  component  b^Q  of  the 
spectral  splitting.  Complete  averaging  is  to  be  expected  at  temperatures  suf'- 
ficiently  high  so  that  — jj— r  «  j,'  whereby  is  the  non-cubic  splitting.  This 
mechanism,  however,  would  cau^e  a  ||mperature  dependent  broadening  of  the  EPR 
lines  to  a  width  of  the  order  of  which  would  mean  that  between  T  =  0 

and  T  such  that  b^Qr/hj  an  appreciable  line  broadening  would  be  observed.  At  the 
latter  temperature  one  would  observe  at  the  most,  one  very  broad  line  centered 
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uzo 

at  the  position  of  the  (-1/2,  l/2)  transition  with  a  width  in  field  AH~  ST  No 

increase  in  the  line  widths  and  no  corresponding  loss  in  intensity  is  observed  in 

++  + 

SrTiOy  Thus,  such  a  random  rotation  is  not  present.  All  the  Gd  ions  (and 
thus  the  whole  lattice)  move  coherently  from  positions  of  tetragonal  to  positions 
of  cubic  symmetry.  This  can  only  be  a  consequence  of  a  cooperative  phenomenon 
due  to  whichllarge:  regions  (containing  many  unit  cells)  of  the  crystal  change 
structure  in  an  ordered  fashion. 

B.  In  this, as  in  other  similar  situations, we  are  faced  with  the  problem  of 

deciding  in  which  site  of  the  host  lattice  the  impurity  substitution  takes  place. 

For  both  Fe+++  and  Gd+++  there  is  a  valence  difference  of  one  unit  to  either 

metallic  constituent  of  the  host.  From  the  EPR  data  in  the  cubic  phase,  we  deduce 

for  both  impurities  that  the  mechanism  responsible  for  the  charge  compensation 

does  not  distort  the  cubic  site  symmetry.  In  the  Gd  spectrum,  a  set  of  low 

intensity  lines  was  observed  indicating  that, for  about  2%>of  the  Gd  ions  in 

these  dilute  crystals,  there  are,  deviations  from  cubic  symmetry  even  at  room 

7 

temperature.  It  was  suggested  that  the  lack  of  cubic  symmetry  in  such  a  spectrum 

was  due  to  association  of  the  Gd  ions  with  A1  impurities.  Analyses  of  these 

crystals  indicate  the  A1  concentration  to  be  far  too  low  to  explain  the  intensity 

of  these  lines.  It  seems  more  probable  that  one  is  dealing  with  charge  compen- 

23 

sation  by  0  vacancies. 

IV 

The  Mn  spectrum,  which  is  the  one  that  should  not  exhibit  the  effects  of 

_| — | — I — f- 

any  charge  compensation  if  we  assume  that  it  substitutes  for  Ti  ,  does  have  an 
axial  field  splitting  in  cubic  SrTiO^.  The  reasons  for  this  are  not  completely 
understood.  We  shall  base  further  discussions  on  the  assumption  that  the  size 
consideration  determines  in  what  site  the  substitution  takes  place.  .Comparison 
of  ionic  radii2^’  25  indicates  that  Gd+++  substitutes  for  the  divalent  alkaline- 
earth,  whereas  Fe  substitutes  for  Ti 

C.  The  temperature  dependence  of  the  Gd  spectrum  indicates  that 
of  the  parameters  in  the  spin  Hamiltonian  only  b£Q  chan'ges  appreciably  with 
with  temperature  below  the  phase  transition.  This  suggests  that  the  mechanism 


12- 


responsible  for  the  b2Q  term  is  independent  of  the  mechanism  responsible  for  the 

other  terms.  Furthermore,  this  mechanism  is  only  present  in  the  tetragonal 

phase.  The  expansion  of  the  crystalline  potential  in  spherical  harmonics  around 

2 

the  alkaline  earth  site  (or,  also,  the  Ti  site)  has  no  V2Q  =  A2Qr  Y(9,  <j>)  (axial) 
component  in  the  cubic  phase.  ^  Such  a  component  will  appear,  however,  with 
the  tetragonal  distortion.  Let  us  assume  that  we  have  a  purely  ionic  crystal. 

Let  us  aj.so  assume  that  the  tetragonal  distortion  is  such  that  the  various  in¬ 
equivalent  sublattices  move  rigidly  in  the  z  direction  and  contract  uniformly  in 
the  transverse  plane.  These  are  the  only  possible  displacements  that  will 

conserve  one  molecule  per  unit  cell  of  the  crystal.  Diffraction  studies  in  BaTiCL 

27  ^ 

indicate  this  to  be  the  case  in  the  tetragonal  phase.  Following  Kenzig  we  call 

these  various  tetragonal  displacements  5z,p.,  6Zq  ,  6Zq  and  6a  the  change  in  the 

transverse  dimension.  Then,  a  point  charge  calc&latioff  shows  that  at  the  alkaline 
♦ 

earth  site 


6V 


20 


”  ,r<r2>4. 


6a  6a  ^  ^zTi 

^-27.27(^)  +  20.  5  (— )  +  .  68 (-^-) 


/nr 


a  I  6zn 

'-3.18  (— 


6z. 

-6.  36 


(4) 


2 

where  r  is  the  electronic  radius  vector  and  the  average  <r  >  has  to  be  taken  over 
the  unfilled  shell  radial  electron  wave  function  (4f  or  3d)  of  the  magnetic  ions,  tr 
is  a  shielding  parameter  which  should  account  for  the  screening  of  the  potential 
by  outer  electron  shells,  e  is  the  electronic  charge  and  a  is  the  edge  of  the  initial 
cubic  unit  cell  in  the  undistorted  phase.  Various  perturbation  schemes  that  per¬ 
mit  calculation  of  the  b.  parameters  for  the  ground  state  of  Gd+++  have  been 
18  *  >  in 

discussed.  They  involve,  in  addition  to  the  crystalline  field,  the  spin-orbit 
or  the  spin- spin. interactions  of  the  4f  shell.  Such  mechanisms  may  involve 
<5^20  either  quadratically  or  linearly.  If  the  important  mechanism  were  quadratic 
in  ^V2(),  one  would  expect  strong  temperature  dependence  to  be  exhibited  also  by 
the  b.Q  parameter,  since  we  have  for  the  square  of  the  second  order  spherical 
harmonic  the  Clebach  Gordan  expansion. 


^20^20  ^(2,  2,  4,  0,  0)Y^0 

i  =  0 


Here  we  have  neglected  the  contribution  to  this  potential  due  to  the  induced 
polarization  on  the  0-"  ions.  This  is  justifiable  in  view  of  the  qualitative 
nature  of  the  argument. 
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Thus  the  assumption  of  a  dominant  mechanism  for  b^Q 
reasonable. 


linear  in  <$V20 


seems 


The  lowest  second  order  perturbation  scheme  linear  in  b2[)  involves 
spin-spin  interactions  whose  magnitude  is  very  difficult  to  estimate. 


The  next  higher  order  in  which  we  have  a  term  linear  in  6V20  is  the 
fourth,  with  a  scheme  that  involves  off  diagonal  matrix  elements  of  the  spin- 
orbit  coupling  and  of  SV20  between  states  within  the  ground  configuration  4f7. 
This  contribution  to  b2^  may  be  written  as 


<8s7/2|l-s)6p7/2><6p7/2|>il-  s|-6d7/2> 


20  r 


E(8s7/2>  -  E(6P7/2I 


e^p7/2)  -  E(6D7/2) 


X 


■f  /2  l«v20 1 6p7/2>  <6p7/2 1 XL'  sl  £’S7/2>> 

1 8 

An  estimate  for  this  effect  has  been  made  for  GdCl3-  For  want  of 
better  information,  we  take  over  this  estimate  from  the  other  crystal  and 
determine  the  ratio 


A20J<r  > 
b20 


-0.  4 


(5) 


Using  displacement  values  obtained  by  diffraction,  ^  we  combine  (4)  and  (5) 
and  under  the  assumption  that  this  is  the  most  important  mechanism  operative 

2  -  ft  # 

in  b2Q,  we  estimate  (<r< r  >)2~.  5x10  cm  which  is  reasonable. 


Unpublished  Fermi  Thomas  calculations  by  G.  Koster  and  H.  Statz  yield,  for 
<r  >2  =  .  53A.  Naturally  the  value  of  the  comparison  is  only  to  show  that  the 
estimated  value  for  Gd+++  has  the  right  order  of  magnitude  and  thus  our 
assumption  of  the  mechanism  for  b2g  is  not  unreasonable,  on  this  count  at  least. 
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f 

J 


Note  that  in  (4)  <5a<0  for  a  contraction.  This  is  really  the  only 
important  term  for  such  a  crude  estimate.  Since  in  BaTiC>3  the  polarized 
phase  may  be  considered  as  resulting  from  a  small  distortion  of  a  cubic 

5  9. 

structure  with  a  center  of  inversion,  the  piezoelectric  strain  —  should 

cl 

depend  quadratically  on  the  permanent  moment  P.  For  small  displacements 

5  z  ^ 

one  may  neglect  the  (  — )  terms  in  (4)  and  thus  we  might  expect,  if  the 

a  2 
assumption  of  the  linear  mechanism  is  correct,  that  b2Q  ccP  .  Incidentally, 

for  rigid  sublattice  displacements  we  see  in  (4)  that  A^q  depends  quadratically 

on  ,5  z.  .  This  is  due  to  the  fact  that  6‘V20  =  0  in  the  cubic  structure  is  an 

extremum  for  such  rigid  sublattice  displacements.  Since  the  permanent 

moment  is  parallel  to  z,  it  should  depend  linearly  on  fi^.ahd  thps,  ev,en  . 

2  c 

with  second  order  terms  included, b^Q  oo  P  .  The  proportionality  b^oc  ^ 
is  indeed  verified,  at  least  with  the  same  degree  of  accuracy  as  that  of 

oc  P2,  as  is  seen  in'F-ig.  7. 28  Here  a  =  V1/3-  where  V  is  the  volume  of 
the  unit  cell.  This  fact  seems  to  support  the  assumption  of  a  mechanism  for 
b2Q  linear  in  6V2q. 


<5V 


D.  Let  us  now  turn  to  the  SrTiO^  data. 


From 


20 


>20 


estimated  above, 


from  the  lattice  constant  a  =  3.  9A  and  from  b2Q  at  4°K  we  estimate  for  a  similar 
mechanism  6a  =  .  01 9A,  whereas  at  77°K  it  would  be  5a  =  .  01 2A.  No  evidence 
for  such  a  large  strain  is  found  in  this  case?^  There  is  an  interesting  point 

J.JLJ. 

however.  In  Fig.  6  we  show  the  temperature  dependence  of  b2Q  for  Gd  in 

this  crystal.  There  is  no  measurable  discontinuity  in  b2Q  or,  for  that  matter,  in 

the  value  of  the  other  parameters  at  the  transition.  Such  would  be  the  situation 

for  a  transition  of  order  n>J.  How  smooth  a  transition  looks  deptends  to  some 

extent  on  the  range  of  the  microscopic  interactions  involved  in  the  phenomenon 

under  observation..  Other  properties  involving  interactions  of  range  shorter  than 

the  electrostatic  potential  might  conceivably  exhibit  much  faster  changes  than-$ 

b2Q.  It  is  a  straightforward  consequence  of  the  "local  field  catastrophe"  theory 

of  transitions  between  a  polarized  and  an  unpolarized  phase  that  if  n  =  2,  the 

dependence  of  the  permanent  moment  on  temperature  below  the  transition  is 
—  30  31 

P  oo  (Ttr  -  T)2.  ’  This  is  true  at  temperatures  sufficiently  high  so  that 

quantum  effects  may  be  neglected.  The  data  in  Fig.  7  seem  to  fit  down  to  50°K 
such  a  law  quite  well  with  T^  =  107.  5°K.  This,  however,  according  to  our 
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assumed  mechanism  for  b20>,  contrary  to  what  happens  in  BaTi03,  would  mean  a 
proportionality  between  b2Q  and  P.  As  a  consequence. linear  terms- in  6z  should 
be  present  in  (4),  which  will  only  be  the  case  if  it  is  assumed  that  the  various 
inequivalent  ,  sublattices  do  not  displace  rigidly  in  the  z  direction  when -the  cubic 
phase  distorts  into  the  tetragonal  phase.  This  in  turn  would  mean  that  in  the 
tetragonal  phase  we  have  a  larger  than  minimum  unit  cell.  The  implication 
would  then  be  that  we  are  dealing  in  SrTiO^  with  a  structural  phase  change  of 
different  nature  than  that  in  BaTiO^  at  120°C  and, since  we  do  have  a  high  and 
T  dependent  dielectric  constant,  the  conjecture  of  an  antiferroelectric  phase 
might  find  partial  justification. 

Using  the  ratio  of  b2^  in  BaTiC>3  at  300°K  to  b2Q  in  SrTi03,  we  estimate 
for  the  latter,  provided  there  is  no  accidental  cancellation in(4),  ionic  dis¬ 
placements  at  4°K  of  the  order  of  5z~  .  0014A  and  at  77°K  6z-»  .  0006A.  They 
would  be  quite  small  at  77°K.  This  would  explain  the  difficulty  in  detecting 
this  phase  transition  by  X-ray  diffraction  methods.  For  this  mechanism  in 
SrTi03,  fc>20<0  indicates  6V2q>0. 

It  is  of  interest  to  compare  the  relative  magnitudes  of  b2Q  in  these 

1J.1  _LJ__L  XJ—f_ 

crystals  for  both  Fe  and  Gd  (the  values  of  these  constants  for  Fe 

jin  BaTiO,  are  less  accurate  but  still  seem  to  exhibit  with  reasonable  accuracy 
2  J 

b20co  P  )  as  seen  jn  Fig.  7. 

b2o(BaTi03)  93Q  l 

Fe+++  -  = -  k 

b20(SrTiO3)  1 6.  1 

b20<BaTlO3>  293.6 

Gd+++  -  =  -  -i 

b20(SrTiO3)  362.  5 

The  temperatures  for  comparison  were  chosen  sufficiently  below  the  transition 
so  that  9b2Q/8T  is  small  and  does  not  affect  the  order-of-magnitude  comparison. 
Similar  ratios  for  the  b^m,  i>  2  in  corresponding  phases  for  both  ions  are  of 
the  order  of  1  independently  of  temperature.  This  indicates  that,  whereas  Gd+++ 
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in  both,  crystals  sees  about  the  same  amount  of  tetragonal  distortion,  such 
distortion  at  the  Fe+++  site  is  much  larger  in  the  barium  salt  than  in. the 
strontium  salt.  If  the  site  assignment  assumed  earlier  is  indeed  correct 
this  would  mean  that  the  difference  in  behavior  at  the  cubic-tetragonal  phase 
transition  of  the  crystals  is  associated  with  a  radical  difference  in  the 
distortion  of  the  Ti-O/ -oetaistedrs?.,  . 

F.  There  is  a  fundamental  problem  with  the  explanation  of  the 

dielectric  behavior  of  SrTiO^  brought  about  by  the  accurate  definition  of 

the  cubic -tetragonal  phase  transition.  Whereas  in  BaTi03  the  actual  phase 

transition  temperature  T^aside  from  hysteresis  effects)  is  very  close  to 

the  Curie  temperature  T  (  «  oc  for  T  >  T  ),  the  two  temperatures 

c  i-ic  32  c 

are  quite  different  in  SrTiO^:.  T£~  30°K,  T^  =  110°K.  The  present  results 

do  not  clarify  this  point.  Nevertheless,there  is  one  aspect  of  these  results 

which  tempts  one  to  some  conjectures  in  this  direction.  If  these  phase 

33 

transitions  are  examined  from  a  lattice  dynamics  point  of  view,  one  may 
describe  them  by  a  change  in  the  relative  equilibrium  position  of  the  ions 
around  which  one  should  expand  the  lattice  potential  in  order  to  satisfy  the 
conditions  for  stable  equilibrium.  As  the  transition  temperature  is  approxi¬ 
mated,  say  in  the  cubic  phase,  the  ionic  vibration  amplitudes  increase  to  a 
point  where  terms  in  the  Hamiltonian  of  order  higher  than  two  in  the  deviations 
from  equilibrium  become  important.  The  ionic  potential  which  had  a  minimum 
for  the  configuration  of  cubic  symmetry  now  developes  another  minimum  for 
the  slightly  different  configuration  of  the  ions,  corresponding  to  a  tetragonal 
unit  cell.  The  transition  temperature  is  the  one  at  which  the  free  energies 

calculated  by  expanding  the  Hamiltonian  around  both  configurations 

33 

become  equal.  It  is  shown  by  Chorcran  that  transitions  such  as  we  have  studied 
may  be  achieved  by  assuming  that  the  modes  in  a  certain  small  range  of  some 
optical  branch  attain  large  amplitudes  and  their  frequency  decreases  to  an 
abnormally  low  value.  In  particular,  such  a  mode  for  BaTiC>3  cbuld  be  an 
optical  mode  with  vibration  along  the  tetragonal  axis  and  in  ,ithe  :  .  ; 

neighborhood  of  wave  vectors  k  =  0  in  the  branch.  The  temperature  dependence 
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of  the  frequency  of  such  a  mode  is  given  by 

w2=K1(T-Tc).  (6) 

This  is  a  consequence  of  the  fact  that  when  the  vibration  amplitudes  are 

sufficiently  large  so  that  higher  than  second  order  terms  in  the  energy 

become  significant,  a  first  order  correction  to  account  for  this  is  the 

assumption  of  spring  constants  K  that  depend  linearly  on  T  (as  in  the 

theory  of  thermal  expansion  )  K  =  Kq+Kj/T).  For  a  ferroelectric  transition, 

K  is  more  than  cancelled  out  by  strong,  long  range  electrostatic  interaction. 

2 

Since  w  oc  K,  a  relationship  of  form  (6)  is  then  obtained.  The  actual  phase 
transition  takes  place  at  a  temperature  slightly  above  Tc  and  thus  u  never 
actually  becomes  zero.  Here  is  where  our  results  might  have  some  signifi¬ 
cance.  If  such  a  lattice  vibration  frequency  would  come  at  all  close  to  the 

frequency  corresponding  to  the  EPR  splitting 

E 

ing  of  the  lines  to  a  width  of  the  order  of  (y) 

maximum  tetragonal  splitting.  No  temperature  dependent  line  broadening  is 

observed  in  SrTi03.  This  would  put  >  103  kMc/sec  at  all  times  which  in 

turn  indicates  a  sizeable  minimum  difference  T,  -T  .  In  the  case  of  BaTiO, 

tr  c  5 

th^s  difference  is  very  small  (~2°K).  The  loss  of  resonance  intensity  observed 

around  the  transition  and  discussed  in  Section  II  might  be  ascribed  to  this  fact. 

There  is  about  a  2°K  hysteresis  associated  with  this  transition  which  means 

that  within  this  2°  region  the  two  phases  exist  simultaneously.  Tnis  might  also 

mean  that  in  the  part  of  the  crystal  still  in  the  cubic  phase  very  large  line- 

widths  could  prevail,  so  large  that  their  contribution  to  the  resonance  lines 

would  disappear  as  id  becomes  sufficiently  low  or  T  comes  sufficiently  close 

to  T  .  If  we  assume  that  the  lines  do  reach  their  maximum  width,  id  could  be 
c 

as  low  as  20  kMc/sec.  However,  much  smaller  broadenings  might  cause 
disappearance  of  the  lines  and  thus  u  minimum  could  be  larger  than  this  by 
at  least  an  order  of  magnitude. 

G.  In  BaTiO,  it  is  assumed  that  the  unstable  modes  correspond  to  wave 

J  2ir 

vectors  k  near  k  =  0.  The  domain  size  is  determined  roughly  by  1  =  y  for 

the  unstable  modes.  In  addition,  modes  where  k  =  0  correspond  to  tetragonal 


it  would  cause  a  broaden- 


—  where  E  is  of  the  order  of  the 

id 
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distortions  with  no  change  of  the  number  of  ions  in  the  unit  cell  (rigid 

distortions  of  the  inequivalent  sublattices).  ■  Distortions  that  correspond  to 

the  increase  in  the  number  of  ions  in  the  unit  cell  by  a  small  factor  would 

come  from  wave  vectors  k  near  the  edge  of  the  Brillouin  zone.  No  domain 

structure  should  be  observed  in  this  case  for  an  unstrained  crystal.  This, 

as  was  said  earlier,  is  a  reasonable  possibility  for  SrTi03.  It  would  mean, 

however,  that  the  size  and  shape  of  the  region  of  uniform  birefringence 

observed  in  SrTi03  would  be  determined  by  the  strain  field  built  into  the 

crystal  at  growth.  This  strain  causes  birefringence  even  in  cubic  phase. 

q 

At  the  phase  transition  the  elastic  constants  change  slightly  and  thus  the 
strain  field  changes.  The  change  in  birefringence  pattern  observed  at  the 
transition  is  a  reflection  of  this  effect. 

V.  CONCLUSION 

The  Gd  EPR  spectrum  has  been  studied  in  detail  in  the  cubic  and 
tetragonal  phase  of  two  isomorphous  crystals  SrTi03  and  BaTi03.  A  bire¬ 
fringence  experiment  is  reported  that  confirms  the  existence  of  such  transition 
as  a  property  of  pure  SrTi03<  This  is  of  some  importance  since  previous 
evidence  of  such  a  transition  came  from  EPR  measurements  of  impurities 
in  the  crystal.  The  temperature  of  the  cubic  to  tetragonal  transition  in 
SrTi03  has  been  determined  to  be  110±2.  5°K.\  This  transition  was  shown 
to  be  of  a  different  nature  than  that  in  BaTiOj.  It  is  a  transition  of  the  order 
of  n>]  and  probably  the  tetragonal  phase  has  no  permanent' polarization,  r  i, 
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TABLE  CAPTION 


Table 


Title 


1 


Spin  Hamiltonian  parameter  of  Gd 


+++ 


in  the  titanates. 


TABLE  1 


Ion 

Host 

T 

g 

b20 

b40 

Gd3+ 

SrTi03 

300°K 

1. 992±. 002 

•  ') 

(-5.  7±.  2)X  lO-4 

Gd3+ 

SrTi03 

77°K 

1. 992±. 002 

(-233.  6±.  5)  X  10-4 

(-4.  8i.  5)  X  10-4 

Gd3+ 

SrTi03 

4.  2°K 

1. 992±. 002 

(-362.  5±.  5) X  10-4 

(-3.  24±.  5) X  1 0-4 

Gd3+ 

BaTi03 

300°K 

1.  995±.  003 

(-293.  6±1 . 0)X10-4 

(+4.  Oil.  0)X10"4 

Gd3+ 

BaTi03 

425<>K 

1.  995±.  003 

- 

~6  X 10-4 

Fe3+ 

BaTi03 

425°K2 

- 

- 

17. 0  X 10-4 

Fe3+ 

BaTiC>3 

300°K2 

- 

b2Q  =  930  X  10"4 

17.  0  X 10“4 

Fe3+ 

SrTi03 

300°K5 

- 

- 

-33  X 10"4 

Fe3+ 

SrTiOs 

77°K5 

- 

b20=7.3X10‘4 

-37  X 10~4 

Fe3+ 

SrTi03 

4.  2°K5 

- 

b2Q=  16.  1  X  10‘4 

-38  X 10-4 

Ion 

Host 

T 

b .  . 

44 

b60 

b64 

C^+ 

SrTi03 

300°K 

- 

(.  5±.  3)X10"4 

- 

Gd3+ 

SrTiOs 

77°K 

(-2.  5±.  3)X10"4 

(-.  25±.  5)X10-4 

t(.  1±.  5)  X  10~4 

Gd3+ 

SrTi03 

4.  2°K 

(-4.  2±.  5)  X 10-4 

(il.  4i-.  5)X10"4 

(-.  69±.  5) X 10  4 

Gd3+ 

BaTi03 

300°K 

(-2.  Oil.  0)X10“4 

(1.  6±1.  0)X10"4 

(~0) 

FIGURE  CAPTIONS 
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Angular  dependence  of  some  of  the  EPR  transitions  for  Gd 
in  SrTiO^,  at  room  temperature  at  a  frequency  of  16  kMc/sec. 
6=0  corresponds  to  H  parallel  to  the  [  100] . 

+++ 

Angular  dependence  of  some  of  the  EPR  transitions  for  Gd 
in  SrTiO^  at  77°K  at  a  frequency  of  16  kMc/sec.  0  =  0 
corresponds  to  H  parallel  to  the  tetragonal  axis. 

Angular  dependence  of  the  EPR  spectrum  of  Gd+++  in  tetragonal 
BaTiO^  at  room  temperature,  with  the  magnetic  field  in  a  (100) 
plane.  The  origin  for  the  angles  is  the  [100]  tetragonal  axis. 
The  frequency  is  1 6  kMc/ sec. 

Temperature  dependence  of  the  Gd+++  EPR  spectrum  in  SrTiO-j, 
with  the  magnetic  field  parallel  to  a  [  100]  direction  and  at  a 
frequency  of  16  kMc/sec.  The  circles  correspond  to  lines  due 
to  the  ions  whose  tetragonal  axis  is  parallel  to  the  fieid. 

Temperature  dependence  of  the  EPR  spectrum  of  Gd+++  ions 
in  BaTi03,  in  the  neighborhood  of  the  tetragonal  to  cubic  phase 
transition,  with  the  magnetic  field  parallel  to  the  tetratonal  axis. 
The  frequency  is  again  16  kMc/sec. 

Temperature  dependence  of  two  of  the  parameters  in  the  spin 
Hamiltonian  for  Gd+++  in  SrTiO^. 

Dependence  on  the  polarization  of  the  piezoelectric  strain  and 
the  second  order  spin  Hamiltonian  parameter  of  Gd+++  in 
tetragonal  BaTi03> 


Solid  lines  -  From  spin  Hamiltonian 
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